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Baryogenesis with QCD Domain Walls 
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We propose a new baryosymmetric mechanism for baryogenesis which takes place at the QCD 
scale and is based on the existence of domain walls separating the metastable vacua from the 
lowest energy vacuum. The walls acquire fractional negative and positive baryon charges, while the 
observed baryon asymmetry is due to a non-zero value of the 9 angle at the temperatures near the 
QCD chiral phase transition. The regions of metastable vacuum bounded by walls carry a negative 
baryon charge and may contribute a significant fraction of the dark matter of the Universe. 
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1 . The origin of the asymmetry between baryons and 
antibaryons, and, more specihcally, the origin of the ob¬ 
served baryon to entropy ratio ns/s ^ 10“® (jib being 
the net baryon number density, and s the entropy den¬ 
sity) is still a mystery. In order to explain this number 
from symmetric initial conditions in the early Universe, 
it is generally assumed that three criteria, first laid down 
by Sakharov must be satihed: (1) there must exist 
baryon number violating processes; (2) these processes 
must involve C and CP violation; and (3) they must 
take place out of thermal equilibrium. 

Many mechanisms for baryogenesis have been proposed 
over the past twenty years (for recent reviews see e.g. 
|^|-Q), but all involve new physics, often at extremely 
high energy scales. A potential problem for all baryo¬ 
genesis mechanisms which operate at scales higher than 
the electroweak scale Aew is that electroweak sphaleron 
effects wash out any pre-existing baryon number unless 
other quantum numbers such as B — L (where B and 
L are baryon and lepton number, respectively) are vi¬ 
olated during baryogenesis. Therefore, in recent years 
electroweak baryogenesis has received a lot of attention, 
see e.g. [H|]. In the standard model, however, CP viola¬ 
tion is too small to generate the observed value of nsjs- 
In addition, given the present bounds on the Higgs mass, 
it appears unlikely that in the standard model the elec¬ 
troweak phase transition is sufficiently strongly hrst or¬ 
der to generate the out-of-equilibrium bubbles required 
to satisfy the third Sakharov criterium. This problem 
can be circumvented within topological defect-mediated 
electroweak baryogenesis but only at the cost of 

introducing new physics above the electroweak scale. 

Recently, it was suggested Q that the Sakharov crite¬ 
ria could be satisfied at the QCD scale. In this paper, we 
propose an explicit mechanism which can explain the ob¬ 
served value of ns/s and which takes place at the QCD 
scale without introduction of any new physics. Our mech¬ 
anism, however, is strictly speaking not a baryogenesis, 
but rather a charge separation mechanism, since the to¬ 
tal baryon current is exactly conserved. In this aspect, 
our approach makes use of ideas proposed in the context 


of toy models in |8||^. 

The first essential ingredient in our scenario is the re¬ 
alization that there exist several vacua in QCD whose 
energy degeneracy is broken only by a very small amount 
AU ~ vciq (rriq is a quark mass) compared to the heights 
of the barriers between the vacua. This leads to the 
domain walls which provide the locus of out-of-thermal 
equilibrium required by the third Sakharov criterion. 

This picture of the vacuum structure is based on 
the analysis of the general properties of the 9 depen¬ 
dence in QCD. It is known that physical values de¬ 
pend on the strong CP parameter 9 through 9/Nf 
(or 9/Nc in pure gluodynamics, Nf and Nc are the 
number of flavors and colors, respectively) and, at the 
same time, they must be 27r periodic functions of 9. 
This is only possible if an effective potential is multi¬ 
valued, and in the partition function one sums over all 
branches [p|Jl^]. This yields the vacuum structure de¬ 
scribed above. In terms of the original theory this pre¬ 
scription means summation over all topological classes 
Z ^ / DAexp [—S'o-I-fQ(0-I-27rn)] imposing the 

topological charge Q quantization. This picture arises 
in supersymmetric and non-supersymmetric the¬ 
ories. The appearance of domain walls in gauge theo¬ 
ries was also advocated by Witten within different frame¬ 
works 0. 

The second crucial input is the fact that the domain 
walls may carry a nonvanishing fermion number, a well- 
known IMl property of many solitonic configurations. 
The domain walls are also locations of maximal CP- 
violation [Q due to the phase difference of the chiral con¬ 
densate across the wall, demonstrating that the second 
Sakharov criterium is satisfied. In what follows we use 
a specific low energy effective Lagrangian for QCD [ p^ , 
though one can expect that the general properties de¬ 
scribed above are inherent in QCD, and should appear 
in any approach. 

Making use of these ingredients in the context of the 
early Universe, it follows by the Kibble-Zurek mechanism 
||l^ , |l6| that at the QCD phase transition at the tempera¬ 
ture Tc ~ 200 MeV a network of domain walls separating 
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different vacua will form. Because of the energy bias 
AE, walls surrounding regions of metastable vacua will 
tend to contract. However, when coupled to fermions, 
they can be stabilized, yielding some compact objects, 
not necessarily of spherical form, to be called B-shells 
in what follows. (This name is due to the fact that the 
(anti-) baryon charge is concentrated on the surface of 
such a compact object. This bears some resemblence to 
the SLAG bag |^.) The compensating positive baryon 
number is left behind in the bulk. The B-shells, in turn, 
can contribute a substantial fraction of the dark matter 
of the Universe. 

2. We start by describing the effective potential 0 
which allows one to analyze the vacuum properties and 
domain walls in QCD. In this approach, the Goldstone 
fields are described by the unitary matrix Uij correspond¬ 
ing to the 75 phases of the chiral condensate: = 

— |(^l'I'/{)| Uij. The effective chiral potential is periodic 
in 9 and takes the form 

Wqcd{Q,U) = - lim — logiy^exp[U£;cos$ 

V^OC V 

1^0 

+ ^VTr{MU + M+U+)]} , (1) 

where V is the 4-volume and 

$ = —-9-\-i- log Det U -\ - 1 , M = d*ag(mi| (^*4'*) |) 

P P P 

{nrii are the quark masses), while E = {bas/{32TT)G‘^) is 
the QCD gluon condensate. Two integers p and q are 
related to a discrete symmetry which is a remnant of the 
anomaly and can be found only by explicit dynamical 
calculations. In the large W limit, q/p ~ l/W in or¬ 
der for the 17(1) problem to be solved. The number of 
metastable vacua is very sensitive to the values of p, q. 
In an approach of 0Jl^ p = 3b = llNc — 2Nf,q = 8 
{Nf is the number of flavors), but here we keep them as 
free parameters. 

One can argue that Eq. i) represents the 

anomalous effective Lagrangian realizing broken confor¬ 
mal and chiral symmetries of QCD. The arguments are 
the following: (a) Eq. (|]) correctly reproduces the Di 
Veccia-Veneziano-Witten effective chiral Lagrangian 0 
in the large W limit; (b) it reproduces the anomalous 
conformal and chiral Ward identities of QCD; (c) it con¬ 
tains the built-in topological charge quantization, which 
shows up in (® via the sum over the integers I and pres¬ 
ence of cusp singularities at certain values of the fields. 

Thus, the effective chiral potential i) satisfies all gen¬ 
eral requirements on the theory and leads to the above 
picture of the vacuum structure and the related domain 
walls. For the domain walls interpolating between the 
true vacuum of lowest energy and the first excited state 
at 0 = 0, the surface energy density a of the wall for the 
potential has been calculated in : 


where = 133 MeV. The energy splitting between the 
ground state and the metastable state at 0 = 0 is 

AE = rUqNf I (^^i) I ^1 - cos + O(m^). (3) 

Due to this energy difference, regions of the metastable 
state tend to disappear, unless they are stabilized by cou¬ 
pling to fermions. In this case, one expects to end up with 
blobs of the metastable state (B-shells). 

3. Next we would like to argue that the above domain 
walls may acquire a fractional baryon charge when cou¬ 
pled to the baryons. As our domain wall configuration 
is a flavor singlet (the q' domain wall [Q), we keep only 
the relevant part of the matrix U = exp[ja(z)], where 
a{z) has a solitonic shape written explicitly in [Q. Ne¬ 
glecting the isospin, we consider the following simplified 
Lagrangian for the nucleon N interacting with the exter¬ 
nal, non-fluctuating chiral field U and additionally via a 
four-fermion interaction: 

£4 = Nidf,-f^N - tun^lUNr - tunNhU+Nl (4) 

—A (NlNr) (NrNl) 

(the choice of the sign of A > 0 corresponds to the repul¬ 
sion in the U(l) channel). 

The important point is that the nucleon mass toat 
takes different values in different vacua: toat = -|- 

mqf{6,k) , where is the nucleon mass in the chi¬ 
ral limit and k labels different vacua. This is related to 
the fact that the chiral condensate varies by an amount 
~ rriq in different vacua. Given this, estimates of a func¬ 
tion f{6, k) are possible but will not be considered here. 
Thus, in the adiabatic approximation to the interaction 
with the domain wall, toat should be considered as a 
slowly varying function of z. Fortunately, its precise form 
is irrelevant for our purposes, only the asymptotia enter 
the final answer, see Eq. i) below. 

Because of the planar symmetry, the problem of com¬ 
puting the charge = f NjqNiUx in our four di¬ 

mensional (4D) theory i) reduces in the mean field 
approximation to the calculation of the corresponding 
charge B^"^^ = f ’ijjjo'ipdz of a two dimensional (2D) the¬ 
ory. Omitting the details of this calculation which will 
be given elsewhere, we here present the final answer for 
the 2D baryon number on the wall: 

z—-\-oo 

, ( 5 ) 

Z— — 00 

where A has dimension of mass and can be found in terms 
of A, and a degeneracy factor n = N/S, see below. 

To find the original 4D baryon charge B^'^\ we should 
take account of a degeneracy related to the symmetry un¬ 
der the shifts along the wall plane. In many body physics 


B^ ' = — arccos 

TT A 
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the definition of the charge is B = J N'^joN^d^x 
where we should sum over all particles (possible quan¬ 
tum states). In our case such a summation yields 


Q = 



dxdydp^dpy ^ n(2)^ 

{2ny 


( 6 ) 


where g = 2 ■ 2 describes the degeneracy in spin and 
isospin. For a fixed number of quantum states N = 
9 ^/ ~ ^ 4 ^^ {S is the area of the wall), the fermi 

energy Ep oi the domain-wall fermions is determined by 


Ef = 



pd?p 
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(7) 


The size of the surface which can accommodate the fixed 
number of fermions N can be found from the minimiza¬ 
tion equation ^^\N=const = 0 ) where the total energy of 
the fermions residing on a surface S is given by 


Eq = aS + 



( 8 ) 


The minimization then relates the density of fermions per 
unit area n = N/S to the tension a (^). We obtain 


Q 


-ST^ai, 


ai = -■ 


j2/3 

w 


(^)l/3^(2) ^ 

'47r'^ 


(9) 


4. Now we are ready to introduce our proposed baryo- 
genesis (charge separation) mechanism. At the QCD 
phase transition at the temperature ~ 200 Mey, 
the chiral condensate forms. Because of the presence 
of nearly degenerate states, a network of domain walls 
will arise immediately after Tc. At 0 = 0, there are two 
degenerate metastable states \B) and \C) above the true 
vacuum of lowest energy |A). For simplicity, we ignore 
metastable states of higher energy, although they might 
be important for evolution of the domain wall network 
[^ . A CP transformation exchanges the states \B) and 
\C). Respectively, the baryon charge of a A — B wall 
will be opposite to that of a A — C wall, and no baryon 
asymmetry can be produced in this case. The situation 
becomes different, however, when 0 7 ^ 0 at the temper¬ 
atures close to Tc- This case will be considered below. 
(It is implied that the strong CP problem will be cured 
by an axion at lower temperatures. At T ^ Tc, the ax- 
ion field does not yet settle in its ground state, and thus 
0{Tc) might be of order unity. Note that as long as an 
initial value 0(Tc) is the same in the entire observable 
Universe, so will be the sign of the baryon asymmetry. 
This occurs if the Universe undergoes inflation either af¬ 
ter or during the Peccei-Quinn symmetry breaking.) In 
this case there is a splitting AE ^ Oruq between the en¬ 
ergies of lowest metastable vacua, which translates into 
the splitting AM ~ 9M for the masses of B-shells with 
negative and positive baryon charges (here M stands for 


the B-shell mass at 0 = 0). If the B-shells reach thermal 
equilibrium at a temperature T, the relative densities of 
negative and positive B-shells will be fixed by the Boltz¬ 
mann formula, and one type of B-shells will thus be very 
strongly suppressed. (A similar interplay between the 
spontaneous and explicit breaking of CB was considered 
in 1211 within a different scheme.) Assuming that this is 
the case, we will therefore discuss the evolution of the 
domain walls of only one (negative) type. 

The initial wall separation depends on the details of the 
chiral phase transition. We will assume that an infinite 
domain wall network will exist until a temperature Td, 
after which it decays into a number of finite clusters of 
the false vacuum. The typical wall separation ^ = ^{Td) 
at this temperature depends on the initial wall separation 
at formation, the details of the damping mechanism and 
the interplay between the energy bias (||) and the surface 


pressure in the walls 


For the time being, we 


will use the value of ^ = ^(Td) as an unknown parameter 
in our scenario, and relegate a brief discussion of possible 
values of ^ to the end of the paper. 

Given the typical wall separation the total area in 
the walls at this temperature is (S') ~ U/^, where 
V stands for the Hubble volume at Td- Note that this 
is parametrically larger by than the Hubble area 

U2/3^ This enhancement means that the baryon charge 
on the wall is in fact not a surface effect, as Eq. (||) 
would naively suggest, but rather a volume effect. Thus, 
the total (anti-) baryon charge stored in the (negative) 
walls at Td is fixed by Eq. (^, while the compensating 
baryon charge is left in the bulk (the positive walls and 
their subsequent decay products are effectively ascribed 
to the bulk). As the entropy density is s = g*T^, where 
g* ~ 10 is the number of spin degrees of freedom in the 
radiation bath at Td, we may estimate the baryon to 
entropy ratio in the bulk at Td'. 


ns 


s Td 


■ 


( 10 ) 


After Td, the domain wall network will break up into 
finite clusters of the false vacuum (“B-shells”). Although 
the question of stability of these B-shells requires a de¬ 
tailed study, qualitatively we expect that the bubbles will 
shrink, but not decay completely since they will eventu¬ 
ally be stabilized by the fermions. As the non-relativistic 
baryons can hardly cross the wall, we expect the shells 
to be stable against the escape of baryons from the inte¬ 
rior, but able to lose heat by baryon pair annihilation and 
emission of the photons and/or neutrinos. The quantum 
stability of the B-shells will be addressed in . Gener¬ 
ally, one expects an exponential suppression of quantum 
decays by the baryon charge of the surface. 

To proceed, we introduce two dimensionless constants 
a 2 and by parametrizing the total area and volume 
of the B-shells as a\VI^ and a^V, respectively. (This 
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parametrization does not imply any specific assumption 
about the form of the B-shells.) In addition, we neglect 
the change of the Hubble radius between the time cor¬ 
responding to Td and the time at which the B-shells are 
formed. We thus obtain 


TlB 2 —1 

— - aia25* 
s 



riB 


s Td 


( 11 ) 


The energy density ps in the blobs of the metastable 
state will redshift as matter, i.e. Pb{T) ~ ■ Hence, 

the B-shells can contribute to the dark matter of the 
Universe. Their contribution to 


PBjteq) _ PB^td) 
Pr{teq) g*TdTeq ’ 


( 12 ) 


(where teg is the time of equal matter and radiation, and 
Pr is the energy density in radiation which, until 
dominates the total energy density) is, roughly. 


Us ~ as 


27r2 


TOgK^'d')! 
n T^T ’ 


(13) 


where we used the expansion of the cosine in Eq. i)- 
This estimate implies that the dominant contribution to 
the energy density pB inside the B-shells is that due to 
the false vacuum energy (^), which is larger than the 
non-vacuum contribution ^ asm^nB with ub given by 
&■ (Here we assume that the baryon density ub is of 
the same order of magnitude in the bulk and interior, 
while the redundant baryons and anti-baryons with zero 
net Ub will eventually annihilate and leave the B-shells 
together with a flux of photons and neutrinos.) 

Comparing ( 0 ) and ( 0 ), we see that the resulting 
baryon asymmetry and contribution of the B-shells to the 
dark matter density are related with each other via the 
geometrical parameters a 2 , as of the distribution of shells 
sizes. To calculate them is a complicated task which is at 
the moment beyond our ability. However, it is of inter¬ 
est to reverse the argument and estimate the parameters 
02 , as assuming that the B-shells do provide the observed 
baryon asymmetry ubIs ~ 10 “® and contribute signifi¬ 
cantly to the dark matter, Us ~ 1 . As can be seen from 
(0) and 0 ), this requires 03 ~ 10 ’’ and ^Td — 10 ®a|, 
if Oi ~ 10 “^ and Td ~ Tc. As 02 < 1 and ^Td ^ 1 
for the thin wall approximation to work, this leaves us 
with the window T~^ ^ for the proposed 

mechanism to be operative. This corridor is consistent 
with the Kibble-Zurek scenario for the topological defects 
formation in the early Universe 0 , ^^ . 

In conclusion, qualitative as our arguments are, they 
suggest that baryogenesis can proceed at the QCD scale, 
and might be tightly connected with the origin of the 
dark matter in the Universe. An alternative scenario 
could also be considered, in which the charge separa¬ 
tion proceeds due to different interactions of particles 


and anti-particles with a CP-breaking domain wall with 
zero baryon charge, similarly to what happens in elec- 
troweak baryogenesis. The blobs of metastable vacua 
might be stabilized in this case by the fermi pressure of 
anti-baryons trapped inside during the evolution of the 
walls. We emphasise that the suggested ideas can be, in 
principle, experimentally tested at RHIC. A more quan¬ 
titative and detailed analysis will be given in 0 ] . 
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